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We study a set of crossed ID systems, which are coupled with each other via tunnelling at 
the crossings. We begin with the simplest case with no electron-electron interactions and find 
that besides the expected level splitting, bound states can emerge. Next, we include an external 
potential and electron-electron interactions, which are treated within the Hartree approximation. 
Then, we write down a formal general solution to the problem, giving additional details for the case 
of a symmetric external potential. Concentrating on the case of a single crossing, we were able to 
explain recent experinents on crossed metallic and semiconducting nanotubes [J. W. Janssen, S. G. 
Lemay, L. P. Kouwenhoven, and C. Dekker, Phys. Rev. B 65, 115423 (2002)], which showed the 
presence of localized states in the region of crossing. 

PACS numbers: 73.21.Hb, 73.22.-f, 73.23.Hk, 73.43.Jn 
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INTRODUCTION 



Physics in ID systems manifests a number of peculiar 
phenomena, such as spin-charge separation, conductance 
quantization,— and anomalous low-temperature behavior 
in the presence of backscattering impurity^ It is reason- 
able to expect that the more complex structures com- 
posed of crossed ID systems, such as crossings and arrays, 
should exhibit some particular features as well. Although 
the transport properties of crossed ID systems and their 
arrays have been thoroughly studied both theoretically 3 
and experimentally^^, the electronic structure of these 
systems is much less understood and the interpretation 
of existing experimental results is challenging. Recent 
scanning tunnelling microscopy (STM) experiments on a 
metallic carbon nanotube crossed with a semiconducting 
one^ have shown the existence of localized states at the 
crossing which are not due to disorder. However, these 
localized states do not appear systematically in all exper- 
iments, i.e. the effect is highly dependent on the nature 
of the carbon nanotubes (metallic or semiconducting), of 
the barrier formed at the crossing, etc. Aiming at clar- 
ifying this problem, we present in this paper a detailed 
study of tunnelling effects between crossed ID systems 
in the presence of potential barriers for massive quasi- 
particle excitations. Because effects of electron-electron 
interactions can be reasonably incorporated in a ran- 
dom phase approximation (RPA),— & we study a simpler 
model, accounting for electron-electron interactions only 
within Hartree approximation. The outline of this pa- 
per is the following: in section II we introduce the model 
that we are going to use to describe the array of crossed 
nanowires. In section III we consider a particular case 
of free electrons and write down explicit solutions for the 
case of one and four crossings. Section IV contains formal 
general solution with additional details given for the case 
of a symmetric external potential. We demonstrate the 
effect of tunnelling on the electronic structure of single 
crossings in Section V and qualitatively discuss different 
possibilities depending on the external potential. Section 



VI contains quantitative analysis and comparison with 
available experimental data of the electronic structure of 
single crossing for different values of parameters. Our 
conclusions and open questions are presented in Section 
VII. 



II. THE MODEL 

We consider a system composed of two layers of crossed 
quantum wires with interlayer coupling. The upper 
layer has a set of parallel horizontal wires described 
by fermionic fields ipj(x), whereas the lower layer con- 
tains only vertical parallel wires described by the fields 
<fii(y). The wires cross at the points (xi,yj), with 
i,j G Z and the distance between layers is d, with 
min(\xi - x i+ i\, \yj - yj+i\) > d, see Fig.l. 




Figure 1: 2D array of crossed wires. 
The partition function of the system reads 

Z = J d[^-]<#*]%*]<%*]e- s / ft , 
with the total action given by 

S = Sq + Ssct + Dint- 



(1) 



(2) 



The first term accounts for the kinetic energy and exter- 
nal potential V° xt (x) , which can be different in each wire 
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and may arise, e.g., due to a lattice deformation, when 
one wire is built on top of another, 

r h/3 



The second term of Eq. (J2j) describes scattering at the 
crossings {xi,yj), 



So = V / dr dxi/j*(x,T)Gj^ J (x,T) 

rh/3 r 

whe 



/■ft/3 
^ Jo 
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Notice that the matrix element Uij describing intra-layer 
contact scattering can, in principle, be different from C7y, 
but both must be real. On the other hand, the con- 
tact tunnelling (inter-layer) coefficient between the two 
Here, fi xy denotes the chemical potential in the upper crossed wires can be a complex number, since the only 



(4) 



(fi x ) or lower (fi y ) layer. 



constraint is that the matrix above must be Hermitian. 



The third term in Eq. ([2]) accounts for electron-electron interactions, 



fitt 



dT 



dr' dx dx'^*(x 1 T)^*(x' 1 T')V c - c {x- x>)tl)j(x,T)ipj{x' ,-r 1 ) 



It. 



fit) 



dr 



fiu) 



dr' / dy / d 1 /tf(y,T)tf(y'S)V e ^(y-rt<pi(y,T)<p i ( V 'y). 



(6) 



III. FREE ELECTRONS CASE 

We start by considering a very simplified case, 
namely, free electrons (no electron-electron interaction, 
V Q ~ c (x) = and no external potential, V° yit {x) = 0). 

Moreover, we assume Uji — Uji — and put [i x = /j, y = 
fi. The interlayer tunnelling is assumed to be equal at 
each crossing point = T and to have a real and posi- 
tive value. In such a case, the partition function consists 
of only Gaussian integrals. We can then integrate out the 
quantum fluctuations, which reduces the problem to just 
solving the equations of motion. Considering a real time 
evolution and performing a Fourier transformation in the 
time variable, we are left with the following equations of 
motion for the fields: 

where m denotes the electron mass and E is the energy 
of an electron state. Firstly, we evaluate the solutions 
for the case of free electrons without tunnelling and then 
we investigate how the addition of tunnelling changes the 
results. The solution for the free electron case consists of 



symmetric and antisymmetric normalized modes, 

ips(x) = —j= cos(k s x), ip a {x) = —j= sin(/c a x), (8) 

respectively. The corresponding momenta k s and k a de- 
pend on the boundary conditions: with open boundary 
conditions k s — ir(2n + l)/2L, k a — irn/L and with pe- 
riodic boundary conditions k s — k a — im/ L for a wire of 
length 2L and n integer. To find the solution for the case 
with tunnelling T ^ 0, we have to solve Eqs. ([7]). These 
equations are linear, therefore, the solution consists of a 
homogeneous and an inhomogeneous parts, 

i> j (x)=rf™(x)+*l>Y h (x), (9) 

which are 

rf om (x) = A je lkx + B je -' lkx , (10) 
l 

Imposing open boundary conditions, ijjj{±L) = 0, we 
find 

A 3 e lkL +B je - tkL + i/jf h {L) = 0, 
A je ~ lkL +B jG lkL +ipf h (-L) = 0. (12) 
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Writing the above equations in a matrix notation and 
inverting yields 



A* 



1 



Bj J 2isin(2£T) \ -e 



ikL g—ikL 
n — ikL „ikL 



Substituting explicitly the expression for ipj(±L) given 
by Eq. (fTTj) and using the mathematical identity 

(Akx p -ikx\( e ikL - e - lkL \ ( sin(£T- kxi) 
v ' ' \ - e - lkL e lkL J \ sin(kL + kxi) 

= cos (2kL) cos(kx — kxi) — cos(kx + kx{) 
leads, after simplifications, to the solution 



i(y) = -T^G{y,yi)il>i(xi), 



where, for open boundary conditions, 



(13) 



G Q (xi, Xj, E) 



[cos(fca;i 4- kxj) 



ft 2 fc sin(2fcL) 1 
cos(2kL - k\xi - Xj\)), (14) 



and the energy E is related to k as E = h 2 k 2 /2m. Similar 
calculations can be performed for the case of periodic 
boundary conditions, yielding Eq. (fl3|) with 



Gp (xi , Xj , E ) 



h 2 ksm(kL) 



cos(kL — k\xi — Xj\). (15) 



A. Two crossed wires 

In particular, for the simplest case of a single horizon- 
tal and a single vertical wires, with just one crossing at 
(xo,yo), the solution is: 



<p{v) 



-TG(x,x ,E)<p(y ) 
-TG(y,ya,E)1>(xo). 



(16) 




Figure 2: Two crossed wires. 

By substituting (x,y) — (xo,yo), we find that at the 
crossing point 



ip(x ) = -TG(x ,x ,E)(p(y ) 
f(yo) = -TG(y ,y ,E)i/j(x ). 



(17) 



The consistency condition requires that 

1 TG{x ai x ,E) 
TG(y ,y ,E) 1 



0. 



(18) 



or 



T 2 G(x , xq, E)G(y , y , E) = 1. (19) 

The solution is even simpler if (x ,y ) = (0,0). Then, 
for open boundary conditions, the symmetric modes are 

ip(0)Tm 



ip(x) 

<p{y) 



h 2 k cos(fcL) 
^(0)Tm 



sin(fcL — k\x\), 
sin(fcL — k\y\), 



H 2 k cos(fcL) 

and the antisymmetric modes are left unchanged in com- 



parison with Eqs. 



Also, 



G(0,0,E) 



mtan(fcL) 
Wk : 



and the secular equation (fl9|) becomes 



TmtSLii(kL) 



I 2 



h 2 k 



= l, 



(20) 



(21) 



which splits into two transcendental equations 

k+ 



Tm < j i . 

-— r ta,n(k^L), 
nr 



_ Tm _ 
k = —p- tan(/c L). 

The first one describes the shifted values of scattering 
states energies, whereas the second equation has an ad- 
ditional bound state solution with E < 0, if T > To = 
H 2 /mL. The appearance of the bound state is exclusively 
due to the presence of tunnelling. For an electron in a 
wire of length 2L = 10 3 nm the corresponding value is 
To = 7.62 x 10~ 5 eV-nm and for quasiparticles the value 
of T) is typically larger, inversely proportional to their 
effective mass. Defining then k = —ik~ and taking the 
thermodynamic limit L — * oo, we find |«| = Tm/h 2 with 
the corresponding bound state energy 



E 



T 2 m 



and the wave function given by 



ip(x) 



(22) 



(23) 



The factor 1/2 instead of l/y/2 comes from the fact that 
now an electron can tunnel into the other wire, where its 
wavefunction (p(0) = -ip(0). Eqs. |22j| and {23]) hold for 
both open and periodic boundary conditions. Since the 
threshold value Tq is quite small, the bound state should 
exist for a typical crossing with relatively good contact. 
However, the energy of the state is extremely small, E ~ 
10~ 8 eV if T ~ To. Qualitatively similar results were 
found by numerical computatio n 10 ' 11 of the ground-state 
energy of an electron trapped at the intersection of a 
cross formed by two quantum wires of finite width. 
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Four crossed wires 



and 



For the case of two wires in the upper and two in the 
lower layers, there are four crossings. In this case, the 
self consistent equations read 





" ipi(xx) ' 




' fi(yi) ' 






ipi(x 2 ) 

■02 (El) 


= M( Xl ,x 2 ,E) 




(24) 




. ^2(^2) . 




. 952(2/2) . 


l 



where 



M(x u x 2 ,E) = -T 



<Pl(V2) 

^2(2/2) 



M{y llV2 ,E) 



( G{x x ,xi,E) G(x 1 ,x 2 ,E) 

G{x 1 ,x 2 ,E) G(x 2 ,x 2 ,E) 

G(xi,x x ,E) G(x 1 ,x 2 ,E) 

G{x 1 ,x 2 ,E) G(x 2 ,x 2 ,E) 



ipi(xi) 
ipi(x 2 ) 

1p2(xi) 

i>i (x 2 ) 



r 



(25) 



(26) 



The secular equation then has the form 

det[M(x u x 2 , E)M{y u y 2 , E)-I] = 0, (27) 

which yields a rather complicated transcendental equa- 
tion (/ is the identity matrix) . The spectral equation for 
bound states E < can be significantly simplified in the 
thermodynamic limit L — > 00. Then, with k = in, for 
both open and periodic boundary conditions, the matrix 
elements become 



G(xi, Xj,E) 



h 2 \ K \ 



-\K(Xi—Xj ) I 



(28) 



and the secular equation in Eq. ((27) has 4 solutions with 
negative energy described by 



E 
E 
E 
E 



Here, a\ 



T 2 m 

T 2 m 

T 2 m 
~ 2K 2 
T 2 m 

\k(X2— Xl)\ 



(1 — ai — a 2 + aia 2 ), 
(1 + ai — a 2 — aia 2 ), 
(1 — ai + a 2 — a\a 2 ), 
(1 + ai + a 2 + aia 2 ). 



«2 



-\<V2-vi)\ an d E 



—H k /2m (notice the implicit dependence of a± and 
a 2 on E). The value of a« depends exponentially on 
the distance between the crossing points. In the limit 
\x 2 — x\\, \y 2 ~ y\\ ~ * 00 the value of ai,a 2 — > 0, which 
correspond to four independent crossings with the bound 
state energy E = — T 2 m/2h 2 , the same value as we found 
in the previous case (see Eq. 1(22)1 ). 



C. A regular lattice of crossed wires 

Consider now a regular square lattice, with lattice 
constant a. Then, one has xi = al and yj = aj. 



From symmetry arguments, the wave functions should 

be ip 3 (x) = ipo(x)e iKyaj and (fi(y) = ipo(y)e lK:cal . After 



substituting them into Eq. ((13)) and using Eq. <|2S[) we 
find 



l/jj(x) 



T( Pa(y 3 )- 



JlK^ l x a 



H 2 K 



iK x a 



sinh(/i.T — Kal x )e 
cosh(«a) — cos(K x a) 



<Pi(y) = - Tip (xi) 



sinh(Kx — n{l x + l)a) 
cosh(Ka) — cos^^a) 

sinh(ftx — Kal y )e lKy 



h 2 n 



cosh(fta) — cos(K y a) 



sinh(Ky — k(I v + l)a) 
cosh(Ka) — cos(Kya) 



where l x , l y £ Z, such that al x < x < a(l x + 1) and al y < 
y < a{l y + 1). Therefore, i>j(x t ) = t/j o (0)e^ K " al+K y a ^ 
and <pi(yj) = <p (Q)e i ( K * ai+K » a fi , with ip o (0) and <p {0) 
related by 



^o(O) = -T 

MO) = -t 



sinh(Ka) 



H 2 k cosh(Ka) — cos(A"a;a) 
m sinh(Ka) 



h 2 K cosh(Ka) — cos(Kya) 
Thus, the spectral equation reads 

{mTf sinh 2 (Ka) 



W)(0), 
MO)- (29) 



1 = 



{h 2 K) 2 [cosh(fta) — cos(K x a)] [cosh(Ka) — cos(A"j,a)] 



By performing an analytic continuation k = ik in Eq. 
((29| . we find an equations similar to the one obtained 
previously by Kazymyrenko and Douqoh^ when study- 
ing scattering states in a lattice. The spectral equation 
describes a band formed by bound states with energies 
—T/a < E < 0. The momenta K x and K y run in the 
interval — 7r < K x a,K y a < tt if T > Tf = 2h 2 /ma 
or inside the region | sm(K x a/2) sm(K y a/2)\ < T/Tf if 
T < Tf. Similar results were calculated^ estimated^ 
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and measured^ in the context of hybridization between 
vertical and horizontal stripe modes in high-Tc supercon- 
ductors. 



IV. A MORE GENERAL CASE 

Now we consider a more general model, which takes 
into account the presence of an inhomogeneous potential 
V™ t (x) arising from possible lattice deformations, and 
includes electron-electron interactions V c ~ c (x), which 
will be treated at a mean field level, within the Hartree 
approximation V^~ e (x). Each crossing (xi,yj) is consid- 
ered as a scattering point with tunnelling Ty and scatter- 
ing potential Uij . The corresponding equations of motion 
then read 

D jx ipj(x) + J~][Uijil>j(xi) + Tijipi(yj)]S(x - x{) = 0, 
i 

Di y ^i{x) + ^{U a ipi{yi)+T^i{xi)]5{y-yi) = 0, 
i 

where 

/), 



h 2 d 2 



' - 2,,,,/,- • l > ; ' 



with Vj(x) = Vf xt (x) + V£j c {x). This model is solved 
most easily through the Green's function satisfying 



D jxi Gj(xi,x 2 ,E) =5(xx-x 2 ) 



with 



Gj(xi,x 2 ,E) = G* j (x 2 ,Xi,E), 

and the corresponding open boundary conditions, 

Gj(x!,L,E) = 0, Gj (x u -L,E) = 0, 

or the periodic ones 

G j (x 1 ,L,E) = G j (x 1 ,-L,E), 
G/(x 1 ,L,E) = G/{x 1 ,-L,E), 

where the prime denotes the derivative with respect to 
X\. Note that we consider real time Green's function for 
a particular wire (not the whole system). The solution 
to the model is 

^A x ) = -^2[Uij^j(xi) +Ti j ipi(y j )]G j {x,xi,E), 
i 

Vi(tf) = <Pi(Vi)+T?Mxi)]Gi(y,yi,E),(30) 
I 

which we require to be normalized 

]T (| \Mx)\ 2 dx + J \My)\ 2 dy^j = 1. (3i) 



The self consistency condition for the value of the func- 
tions at crossing points {xi,yj) yields the equations 



E) + 5u)ipj(xi) 



+Ti j G j (x i ,xi,E)tpi(y j )] = 0, 
G s (yj,y h E) + 5ji)(fi(yi) 

i 

+T* l G i (y j ,y l ,E)^ j (x i )} = 0. (32) 

To find nontrivial solutions for the fields ipj{x) and <fi(y), 
the system of homogeneous equations in Eq. <f32|) has to 
be linearly dependent and hence the solution is repre- 
sented by the null space of the system. This means that 
after writing the equations in a matrix form, the deter- 
minant of the matrix should be zero, thus leading to a 
spectral equation for E. Moreover, bound state solutions 
in the thermodynamic limit L — > oo satisfy both open 
and periodic boundary conditions, since i/;(±L) — > and 
ip'(±L) -» 0. 

To understand better the dependence of the Green's 
function Gj{xi,xi,E) on E, we represent the function 
through the solutions of the homogenous equations, 



D jx ipj(x) = 0. 



(33) 



We omit the index j in what follows for simplicity. The 
most general and common representation, which holds 
for any static potential, reads as follows: 



G{x\,x 2 ,E) 



n 



E 



(34) 



Here, the function ip e (x) is the solution of the homoge- 
nous equation 



0. 



(35) 



and the spectrum {e n } is obtained by imposing the corre- 
sponding boundary conditions. Notice that in the present 
representation of G(xi,x 2 , E) the functions ip erl (x) have 
to be orthonormal. By writing G(x\, x 2 , E) in the form 
given in Eq. (|34|) . the following identity arises 



J dx'G(x 1 , x',E)G(x',x 2 ,E) = 



dG(x!,x 2 ,E) 
dE 



(36) 



The case x\ — x 2 = for free electrons is illustrated 
in Fig. 02 where Eq. ff20|) is plotted. If some external 
potential is present, G(xq,xq, E) has the same form but 
the positions of the poles are shifted and the correspond- 
ing values are different. If no regularization is used, the 
calculations for E > must be performed in the finite 
size limit, otherwise with L — > oo the energy distance 
between different modes vanishes and the poles situated 
on the real positive half axis merge to form a branch 
cut singularity. This behavior can be readily seen on the 
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G(0,0,E) 



Figure 3: G(0, 0, E) in units of m/h 2 versus E in units of 
h 2 /2mL 2 . 



example of Eq. f20|) , where can perform an analytic con- 
tinuation, considering k — > k + ik' . Then, in the limit 
L — > oo, tan(fcL + ik' V) = isgn(fc'), and the function 
G(xq, xq) changes sign as one goes from the upper to the 
lower complex half plane for k ^ 0. 

Now we represent the Green's function through the 
solutions of the homogenous equation 



2m dx 2 



+ V(x) - E ip(x) = 



(37) 



This is a second order differential equation, therefore, it 
should have two linearly independent solutions, which we 
call ipi(x) and ip2(x). Then the Green's function is 



G(xi,x 2 ,E) 



A-i/>i(xi) + B_ip 2 (xi), xi < x 2 
A+ipi(xi) + B + ip 3 (xi), x\ > x 2 



(38) 

where the expressions for the coefficients A- , B- , A + , B + 
(functions of X2), are derived in the Appendix A. In par- 
ticular, for a symmetric potential V(x), we can choose 
a symmetric ip s (x) and an antisymmetric ip a ( x ) solu- 
tions as linearly independent, i.e., ipt(x) = ijj s {x) and 
ip2(x) = ip a (x). Thus we find 



G(x,0,E) 



S 2 Va'(0) 



1ps(x) 1pa{\x\) 



and 



G(0,0,E) = 



ML) ML) 
mMO) ML) 

h 2 M(o)ML)' 



(39) 



(40) 



To obtain the results in the thermodynamic limit L — > oo, 
it is useful to rewrite G{x\,x%) using quantities which do 
not depend on L explicitly. For example, 



G(x,0,E) = G(0,0,E) 



ip s {x) _ m ip a {\x\) 

Mo) v M(o) ' 



(41) 



After substitution of Eqs. (J8J) into Eq. (|38|) and simplifi- 
cation, for the case of noninteracting electrons we find 



G{x\,X2,E) 



H 2 ksm(2kL) 
cos(2kL — k\xi 



cos(/cxi + kx 2 ) 



which is the same expression as in the previous section 
(see Eq. {Hp . This is a posteriori justification of the use 
of the same letter G(xi,X2,E) in the first section. The 
case of a harmonic potential is considered in Appendix 
B. 



V. A SINGLE CROSSING 

Now we apply our results including tunnelling and ex- 
ternal potential to the simpler case of only two crossed 
wires, aiming to compare our findings with experiments. 
Using the general solution given by Eq. lf30j) , and consid- 
ering T — T*, we can write 

ip{x) - -[Ui>(x )+Tip(y )]G 1 (x,xo,E), 
tp{y) = -[U v(y )+TMo)}G2(y,y ,E). 



By substituting (x, y) 
crossing point 



(xoiUo), we find that at the 



[1 + UGi(x ,x ,E)]ip(xa) + TGi(x ,x ,E)(p(y ) = 0, 
[l + UG2(y ,yo,E)My )+TG2(yo,yo,E)i;(xo) = 0. 



The consistency condition requires that 

1 + UGi(xq,x ,E) TGi(x ,x Q ,E) 
TG 2 (y ,yv,E) l + UG 2 (y ,yo,E) 



= 0, (42) 



= [l + UGi(x ,xo,E)][l + UG2(y 0) yo,E)] 
- T 2 G 1 (xo,x ,E)G 2 (yo,yo,E). 

The meaning of this equation becomes clearer in the 
symmetric case, when U = U and G\{xq,xq,E) = 
G2(yo,yoi E) = G. In this case, it reduces to a quadratic 
equation, which bears two solutions, 



G A 



-1 
U + T' 



G- 



-1 
U-T' 



Notice that they differ by the sign in front of the tun- 
nelling amplitude T, which is shifting the potential U. 
Such symmetry effectively reduces the problem to ID 
with effective potential U e g8(xo). Hence, we have 



ip(x ) = <p(y ), 
ip(x ) = -<p(y ), 



M s = u 



U. 



cff 



T, 
U-T. 



(43) 



The shift of the energy levels in a wire due to the presence 
of the 5 potential can be visualized with the help of the 
Green's function expansion, where one has 



g(x ,x ,e) = 



\MJMf 

e n -E 



TT- (44) 



X2\ 



In the case with U e g = 0, the energies are exactly those 
of the poles and, therefore, remain unshifted. How- 
ever, since G(xq, Xq, E) = —1/U e s, the curve actually 
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describes how the energies of the modes change as we 
keep increasing —1/U c g from — oo if [/ c ff > or decreas- 
ing — 1/U c g from +00 if U c « < 0. In the latter case, we 
can run into the region with E < 0, which would corre- 
spond to the appearance of a bound state. Nevertheless, 
to obtain an exact solution, it is more convenient to work 
with the expression for G(xq, xq, E) in terms of the wave 
functions, 

G(0 '°' £) -/«jWl)-^' (45) 

where we assumed x — for simplicity. 



the Fermi energy, thus suggesting that the electron-hole 
symmetry is broken by the presence of some external po- 
tential. The latter may appear due to lattice distortions 
and the formation of a Schottky barrier at the contact be- 
tween the nanotubes i 16 i 17 In the following, we show that 
if the potential is strong enough, localized states can form 
in the spectrum of the massive mode corresponding to the 
van Hove singularity with energy e = £ v b_ — E. Therefore, 
the observed localized states should have Eq = —0.5 eV 
and E x = -0.2 eV. 

To incorporate in a more complete way the effects of 
the Schottky barrier and lattice deformation, we assume 
V cxt (x) to have a Lorentzian shape, 



VI. COMPARISON WITH EXPERIMENTS 




Figure 4: Voltage versus length diagram, which shows the ex- 
perimentally observed density of states. Notice the existence 
of two localized states in black. (Extracted from Ref. 0). 

Now, we will compare our theoretical findings with ex- 
perimental results. We concentrate mostly on the analy- 
sis of a system consisting of two crossed single wall carbon 
nanotubes (SWNTs): a metallic on top of a semiconduct- 
ing (MS) one.— In its unperturbed state, the band struc- 
ture of a SWNT can be understood by considering the 
electronic structure of graphene. Due to its cylindrical 
shape, the transverse momentum of one particle excita- 
tions in a SWNT has to be quantized, whereas the lon- 
gitudinal momentum may vary continuously. Combining 
this condition with the assumption that the electronic 
structure is not very different from that of graphene, one 
finds two different situations, depending on the topol- 
ogy of the SWNT: there are no gapless modes and the 
nanotube is semiconducting, or two gapless modes are 
present and the nanotube is called metallic. Analyz- 
ing the spectroscopic measurements performed along the 
metallic nanotube (see Fig. [4]) and comparing with the 
unperturbed electronic structure, one notices two main 
changes. First, a small quasi gap opens around the Fermi 
energy level ef between ef — 0.2 eV and ef + 0.3 eV in 
the spectrum of the massless modes (corresponding to 
zero transverse momentum). Second, two peaks are visi- 
ble at £0 = Ef — 0.3 eV and e\ = Ef — 0.6 eV in the region 
around the crossing, corresponding to localized states be- 
tween the Fermi energy and the van Hove singularity at 
£vH = £ f — 0.8 eV. Such states are not visible above 



yext M = -TTxW (46) 

Firstly, we study the influence of this potential alone 
on the electronic structure, i.e. we assume that there 
is no tunnelling T = 0, and no electron-electron inter- 
actions. Exact numerical solution of the Schrodinger 
equation shows that an approximation of the potential 
in Eq. (|46"1) by the harmonic one does not change the so- 
lution qualitatively. Therefore, we consider V cxt (x) w 
—V(l — x 2 /b 2 ), which describes a harmonic oscillator 

with frequency lo = J 2V /mb 2 and corresponding spec- 
tra E n = -V + (n + l/2)^2h 2 V/mb 2 for E n < 0. 
Moreover, it is reasonable to assume that the strength 
of the barrier V is of the same order as the energy of 
the bound states and that the potential is localized on 
the same length scale as the localized states. Hence, we 
take V = 0.7 eV and 6 = 4 nm. It follows then from 
our calculations that the difference between neighboring 
energy levels is quite small and there are many bound 
states present in the case when m is the actual electron 
mass. However, assuming m to be an effective electron 
mass, with m — 0.025 m e , which is of the same order 
as the experimentally estimated values m = 0.037 mj^ 
and m — 0.06 m e ,~ we find exactly two pronounced 
bound states: the first one has E = —0.5 eV and is de- 
scribed by the symmetric wavefunction ips(x) as shown 
in Fig. [H whereas the other has E = —0.2 eV and is 
described by the antisymmetric wavefunction ip a (x), see 
Fig. EE Considering Fig. [5J we observe that the localiza- 
tion size of the state is around 10 nm, which agrees well 
with the experimental data. On the other hand, the state 
shown in Fig. [6] has a zero value exactly at the crossing 
and is rather spread, a behavior which is not observed 
experimentally. Besides these two, a number of other 
states are also present in the vicinity of the van Hove 
singularity with E > —0.1 eV. 

Secondly, we take into account electron-electron inter- 
actions to consider other possibilities to obtain two pro- 
nounced bound states. Unfortunately, our approach only 
allows us to incorporate electron-electron interactions at 
the mean-field level by using the Hartree selfconsistent 
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Figure 6: ip a (x)(nm 1//2 ) versus x (nm). 

approximation 

V H~ e ( x ) = I dx'V c -°(x- x')n(x'), (47) 



where n(x) is the electron density, given by 

n(x) = ^\i}} k {x)\ 2 n F {e k - fi). (48) 
k 

Here the summation k goes over energy levels and np-(e) 
is the Fermi distribution. Although it is known that in 
ID systems quantum fluctuations play an extremely im- 
portant role, we nevertheless start with the mean-field 
approximation as a first step to incorporate them in RPA. 
Moreover, we believe that their presence does not qual- 
itatively change the obtained results. To render the nu- 
merical calculation simpler, we consider a delta-like in- 
teraction potential, which leads to 



v h e ( x ) = Von (a;), 



(49) 



By estimating the effective interaction strength Vq ~ 
2itHvf from the Luttinger liquid theory, we obtain that 
V — 3.4 eV-nm for v F = 8.2 x 10 7 cm/s3 Suppose 
that the lowest energy state with E — —0.5 eV is occu- 
pied by an electron with a certain spin. Then, there is a 
possibility to add to the same state an electron with an 
opposite spin. However, due to the repulsive Coulomb 



interaction the energy of the two-electron state becomes 
E = —0.2 eV for Vq = 3.15 eV-nm. The correspond- 
ing self consistent solution is presented in Fig. [71 The 




Figure 7: ip a (x)(nm ' ) versus x (nm) 

state has the same shape as in Fig. [5j but is a bit more 
spread. By comparing the density of states (DOS) dis- 
tribution with scanning tunnelling spectroscopy (STS) 
data for the crossing,— we observe that the inclusion of 
electron-electron interactions (Fig. [Jj provides a much 
better agreement between theory and experiment for the 
E = — 0.2 eV bound state than in the previous case (Fig. 

EJ). 

Thirdly, we take into account tunnelling between the 
wires. Qualitatively, this leads to the splitting of energy 
levels and redistribution of charge density in the wires, 
thus effectively reducing the strength of electron-electron 
interactions. Since we have no information about the 
electronic structure of the semiconducting nanotube, to 
make a quantitative estimation we assume that the effec- 
tive mass is equal in both wires and that the potential is 
also the same. In such a case, from symmetry arguments 
the electron density should be evenly distributed in both 
wires even for a very weak tunnelling. Therefore, the 
electron-electron interactions should be twice stronger 
than in the case without tunnelling, namely, Vb = 6.3 
eV-nm to achieve the same energy value. Moreover, if 
the tunnelling coefficient is large enough, the splitting 
of the energy levels becomes significant and detectable. 
We can estimate the coefficient T, if we assume that 
it has the same order for SM, metallic-metallic (MM), 
and semiconducting-semiconducting (SS) nanotube junc- 
tions. The SS and MM junctions have Ohmic voltage- 
current dependance, characterized by the conductance G. 
Moreover, we can estimate the transmission coefficient of 
the junction as G/Gq ~ (T/2irhvF) 2 , for G/Gq <C 1. 
For MM junctions experimental measurements^ typi- 
cally yield G/Gq ~ 10~ 2 , thus corresponding to T ~ 0.34 
eV-nm. For example, for T = 0.28 eV-nm and V = 0.44 
eV in Eq. (|46|1 . without electron-electron interactions we 
find that the system has two bound states. The lowest 
energy bound state with E = —0.5 eV is shown in Fig. [HI 
Compared with Fig.[U the state has a peak exactly at the 
crossing, corresponding to a local increase of the DOS. 
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Figure 8: ip-(x)(nm 1//2 ) versus x (nm 



The other bound state with E = —0.2 eV is shown in Fig. 
[H Contrary to the previous case, the state has a deep 
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Figure 9: i|i+(i)(nm 1 ^ 2 ) versus x (nm). 

at the crossing, corresponding to a local decrease of the 
DOS. However, these local change in DOS is too small 
to be observable in the present experimental data. If we 
now include electron-electron interactions with Vq = 3.15 
eV-nm and add a second electron with different spin to 
the system, we find that the new state has E = —0.267 
eV and acquires the shape shown in Fig. [I0l 



w(x) 




if the potential in the metallic SWNT is significantly de- 
creased due to screening effects but a Schottky barrier in 
the semiconducting SWNT can reach considerable values, 
sufficient for the formation of the bound states, then the 
latter are also going to be present in the metallic SWNT 
due to tunnelling between SWNTs. Secondly, there is still 
a possibility to find a bound state existing purely due to 
tunnelling, i.e., without external potential, as was shown 
in Eq. ([23|) . and a second bound state may arise with 
different energy due to Coulumb repulsion between elec- 
trons with different spins. However, this is most proba- 
bly not the case we have in the experiments, because due 
to electron-hole symmetry such states would exist also 
above the Fermi energy, a result which is not observed 
experimentally. 



VII. CONCLUSIONS 

We presented several possibilities to explain the ob- 
served localized states at the crossing of metallic and 
semiconducting nanotubes.— All of them require the ex- 
istence of an external potential in the metallic and/or 
semiconducting SWNT to break the electron-hole sym- 
metry, since the localized states were seen only below 
the Fermi energy. Most probably, such a potential comes 
from a Schottky barrier and the effect of lattice distor- 
tions is minimal, since such localized states were, up to 
now, observed only for MS crossings and not for MM or 
SS ones. Moreover, the effective mass of quasiparticle 
excitations should be of order m = 0.025 m e , where m e 
is the actual electron mass, to generate only a few bound 
states localized on a region of approximately 10 nm with 
energy of order of 0.5 eV. The best agreement with the ex- 
perimental data is obtained by assuming that the second 
bound state has a different energy due to the Coulumb 
repulsion between electrons with different spins. The role 
of tunnelling in the observed electronic structure is not 
clear and allows for many interpretations. To avoid such 
ambiguity, the electronic structure of the semiconducting 
nanotube should be measured as well. Moreover, to be 
sure that the available STS measurements indeed repre- 
sent the electronic structure of the nanotube and are free 
of artifacts introduced by the STM tip2i several measure- 
ments with different tip height should be performed. 
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Appendix A 



The last result suggests that there are yet other pos- 
sible interpretations of the experimental results. Firstly, 



Here we consider the Green's function as a function 
of one variable X\ and fix X2 for a moment. Since 
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G(xi, X2, E) is the Green's function, we require it to sat- 
isfy proper boundary conditions G(±L,X2,E) — 0, be 
continuous G(x2 — 0, X2, E) = G(x 2 + 0, X2, E), and also 
G'(x2-0,x 2 ,E)-G'{x2+0,x 2 ,E) = 2m/h 2 . Substitut- 
ing Eq. I|38p into the above requirements one finds 



P 



where 



P = 
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(Al) 



(A2) 



Multiplying the Eq. (|A1[) by the matrix P 1 we find 



(A, 
A. 



where 



C 



/ Ml)\M-l)M x 2) - M-L)M x 2)} 
M-l)[Ml)M x 2) - ML)MM\ 
-Ml)[M-l)MM - M-l)M x 2)] 
\ -M-l)[ML)MM - ML)M x 2)] 



c 



2m 



\Ml)M-l)-M-l)Ml)]- 1 - 



h 2 Wr 

The Wronskian 

W r = MMMM ~ M X 2)M X 2), 

is nonzero for linearly independent functions and its value 
does not depends on the point #2- 



potentials ip(—x) ls a ^ so a solution and both of them are 
linearly independent. Furthermore, we can then compose 
a symmetric ip s (x) — (ip(x) + tp(—x))/2 and an antisym- 
metric tp a ( x ) — — V-'(~ a; ))/2 solutions. In particu- 
lar, for a harmonic potential V(x) = mw 2 i 2 /2, one can 
find such a solution 



(Bl) 



where H [y, x) is the Hermite polynomial for integer v. It 
follows then that 



(B2) 



and 



^(0) = -2wy 



2'Kiom 



x V 4 fruj) 



Moreover, in the thermodynamic limit L — > oo, 



G(0,0,£) = - 



2»V«ftr($-&) 



(B3) 



(B4) 



The Eq. (|B4|1 approaches asymptotically the expression 
for free fermions, as u> — * for E < 0, 



Appendix B 

Suppose that Eq. f37j) has a solution ip( x ) which is nei- 
ther symmetric nor antisymmetric. Thus, for symmetric 
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